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Abstract A scheme for analyzing Thomson scattering of x-rays by warm dense 
matter, based on the average-atom model, is developed. Emphasis is given to x-ray 
scattering by bound electrons. Contributions to the scattered x-ray spectrum from 
elastic scattering by electrons moving with the ions and from inelastic scattering 
by free and bound electrons are evaluated using parameters (chemical potential, 
O ' average ionic charge, free electron density, bound and continuum wave functions, 

and occupation numbers) taken from the average-atom model. The resulting scheme 
provides a relatively simple diagnostic for use in connection with x-ray scattering 
" J5 . measurements. Applications are given to dense hydrogen, beryllium, aluminum, ti- 

tanium, and tin plasmas. At high momentum transfer, contributions from inelastic 
. scattering by bound electrons are dominant features of the scattered x-ray spectrum 

for aluminum, titanium, and tin. 



1 Introduction 

Measurements of Thomson scattering of x-rays provide information on tempera- 
tures, densities and ionization balance in warm dense matter. Various techniques 
for inferring plasma properties from x-ray scattering measurements have been de- 
veloped over the past decade H-22]. Many of these techniques, together with the 



underlying theory, were reviewed by Glenzer and Redmer in Ref . 12311 . 

The present average-atom scheme is based on a theoretical description of x-ray 
scattering proposed by Gregori J^t], the important difference being that parame- 



W. R. Johnson 

University of Notre Dame, Notre Dame, IN 46556, e-mail: johnson@nd . edu 
J. Nilsen 

Lawrence Livermore National Laboratory, Livermore, CA 94551 e-mail: nilsenl@llnl . gov 

K. T. Cheng 

Lawrence Livermore National Laboratory, Livermore, CA 94551 e-mail: ktcheng@ llnl . gov| 



1 



2 



W. R. Johnson, J. Nilsen and K. T. Cheng 



ters used to evaluate the Thomson-scattering dynamic structure function are taken 
from the average-atom model. The average-atom model used here was introduced 



in Ref. 12411 to study electromagnetic properties of plasmas. The scheme developed 



here to analyze Thomson scattering is closely related to that of Sahoo et al. |25fl. 



where a somewhat different version of the average-atom model was used. Predic 



tions from the present model differ substantially from those in Ref. 02511 . The origin 
and consequences of these differences will be discussed later. 

The Thomson scattering cross section for an incident photon with energy, mo- 
mentum (hcoo, liko), and polarization Co scattering to a state with energy, momentum 
(ha>i , Tiki), an d polarization B\ is 

da ( da\ co\ 



where 
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The dynamic structure function S(k, co) appearing in Eq. (Q~|i depends on two vari- 
ables: k = \k-Q — k\\ and CO = coq — CO\. As shown in the seminal work of Chihara 
co) can be decomposed into three terms: the first Su(k, co) is the con- 



tribution from elastic scattering by electrons that follow the ion motion, the second 
S ee (k, co) is the contribution from scattering by free electrons, and the third St,(k, co) 
is the contribution from bound-free transitions (inelastic scattering by bound elec- 
trons) modulated by the ionic motion. The modulation factor is ignored here when 
evaluating the bound-free contribution. For the bound-free scattering, calculations 
are carried out using both average-atom final states and plane-wave final states. Sub- 
stantial differences are found between average-atom and plane-wave calculations, 
particularly in the low-momentum transfer region of the scattered x-ray spectrum. 

The average-atom model is discussed in Sec. [2] followed by a discussion of the 
three contributions to the structure functions in Sec. [3] In Sec. [4] applications are 
given to hydrogen, beryllium, aluminum, titanium, and tin plasmas. 



2 Average-Atom Model 

The average-atom model is a quantum mechanical version of the temperature- 
dependent Thomas-Fermi model of a plasma developed sixty-three years ago by 



Feynman, Metropolis and Teller 112811 . In this model, the plasma is divided into neu- 
tral Wigner-Seitz (WS) cells (volume per atom V ws =A/ pNa, where A is the atomic 
weight, p is the mass density, and Na is Avogadro's number). Inside each WS cell 
is a nucleus of charge Z and Z electrons. Some of these electrons are in bound states 
and some in continuum states. The continuum density is finite at the cell bound- 
ary and merges into the uniform free-electron density n e = Zf /V ws outside the cell, 
where Zf is the number of free electrons per ion. Each neutral cell can, therefore, be 
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Fig. 1 Main plot: The radial density 47lr 2 n(r) in a.u. for metallic Al at k B T = 5 eV (solid black 
curve) integrates to Z = 13 for r < R ws = 2.991 «o- The continuum contribution ATtr 2 'n c (r) (dashed 
red curve), also in a.u., integrates to 3 for r < S ws . The bound Is, 2s, and 2p shells are completely 
occupied at this temperature. Inset: The dashed red curve illustrates the Friedel oscillations of the 
continuum density and shows how N c (r) = n c (r)V W s converges to Zf = n e V ws (solid green line) for 
r > R vs - The chemical potential predicted by the model is /i = 0.2406 a.u. and the number of free 
electrons per ion is Zf = 2. 146. 



regarded as an ion imbedded in a uniform sea of free electrons. To maintain over- 
all neutrality, it is necessary to introduce a uniform (but inert) positive background 
density Zf/V vs - The model, therefore, describes an isolated (neutral) ion floating in 
a (neutral) "jellium" sea. 



The quantum-mechanical model here, which is discussed in Ref. 12411 . is a non- 
relativistic version of the Inferno model of Liberman 02911 and the more recent Pur- 
gatorio model of Wilson et al. oil : it is similar to the nonrelativistic average-atom 
model described by Blinski and Ishikawa Dill . Specifically, each electron in the ion 
is assumed to satisfy the central-field Schrodinger equation 



P Z z 
I---+V 
2 r 



Va{r) = e a Y a (r), 



(3) 



where a = (n,l) for bound states or (e,l) for continuum states. Atomic units (a.u.) 
where e = h = m = 4-neo = 1 are used here. In particular, 1 a.u. in energy equals 2 
Rydbergs (27.21 1 eV), and 1 a.u. in length equals 1 Bohr radius «o (0.529 A). 
The wave function \j/ a (r) is decomposed in a spherical basis as 



Va(r) = ^P a (r)Y[ ama (?)x<, a , 



(4) 



where F; m (r) is a spherical harmonic and % a is a two-component electron spinor. 
The bound and continuum radial functions P a (r) are normalized as 
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/ drP n i{r)P n ii(r) = 8 nn i, (5) 
Jo 

drP el (r)P en (r) = 8(e-e'), (6) 



respectively. The cent ral p otential V(r) in Eq. (0 is taken to be the self-consistent 
Kohn-Sham potential HI 



V(r) =4-71 1 —r' 2 n(r')dr' - 

r > 



K 



1 



(7) 



where the first term in the right-hand side is the direct screening potential with 
r> = max(r, r 1 ) and the second term is the Kohn-Sham exchange potential. While 
electron-electron interactions inside the Wigner-Seitz cells are reasonably well ac- 
counted for by this simple model, it should be noted that eigenvalues in the Kohn- 
Sham potential are poor approximations to ionization energies, leading to inaccurate 
thresholds and peaks of bound-free contributions to S(k, CO), which can differ from 
experiment by 20 - 30%. The electron density n(r) in Eq. © has contributions from 
bound states «t(r) and from continuum states n c {r), 

n(r)=n b (r)+n c (r). (8) 

The bound-state contribution to the density n^(r) is 

4^ 2 «,(r)=^ 2 ^ + ?,w. Tl W 2 , (9) 

~ l+exp[(£„/-ju)/^r] 

where £„/ is the bound-state energy, jj, is the chemical potential, and the sum over 
(n,/) ranges over all bound subshells. The continuum contribution to the density 
n c (r) is given by 

4nr 2 n c (r) = V Hde 2 j 2/ + 1) w T P £l (r) 2 . (10) 

Finally, the chemical potential jj. is chosen to ensure charge neutrality inside the WS 
cell: 

/ A it, J-,, 



Z= / «(r)t/ j r= / 47Zr z n(r)dr . (11) 

Jr<Rv/ S JO 

Equations OHTTJ above are solved self-consistently to give the chemical potential 
jj., the potential energy function V(r) and the electron density n(r). 

The upper limit in Eq. (O is determined by systematic trial and error. The val- 
ues of n and / are increased starting from n = 1 and / = 0. If a state is bound, it 
is included in the sum, otherwise not. At metallic densities and temperatures below 
100 eV (warm dense matter) fewer than a dozen states typically bind. To carry out 
the sum-integral in Eq. ( fTOb for the continuum density, we typically use 12 partial 
waves (I) and 40 to 50 energy points (e) for each partial wave. The energy grid for 
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Table 1 Aluminum at density 2.70 g/cc and k B T = 5 eV. Bound-state and continuum partial-wave 
occupation numbers inside the WS sphere are given, along with bound-state eigenenergies. The 
Ne-like core is seen to be almost completely occupied. The sum of the bound-state and continuum 
occupation numbers is precisely Z = 13. 





Bound States 


Continuum 


State 


occ# e(eV) 


I 


occ# 


Is 


2.0000 -1485.07 





0.9130 


2s 


2.0000 -92.16 


1 


1.3263 


2p 


5.9998 -54.87 


2 


0.6192 






3 


0.1173 






4 


0.0209 






5 


0.0031 






6 


0.0004 






7 


0.0000 


N b 


9.9998 


N c 


3.0002 



the integral in Eq. ( TTOb is chosen using a modified Gauss-Laguerre scheme. Thus, 
one is faced with solving a system of roughly 500 coupled second-order differen- 
tial equations. These equations are solved iteratively using a predict-correct scheme 
based on Adam's method I33I Chap. 2.3]. 

The boundary conditions used in solving Eq. (O deserve some mention. Bound- 
state wave functions and their derivatives are matched at the boundary r = R vs to 
solutions outside the WS sphere (where V = 0) that vanish exponentially as r —> °°. 
Similarly, continuum functions and their derivatives are matched to phase-shifted 
free-particle wave functions at r = R m , It should be noticed that the continuum den- 
sity n c (r) inside the WS sphere, which oscillates as predicted by Freidel 13411 . is 
distinctly different from the uniform free electron density n e . In the present model, 
n c (r) smoothly approaches n e outside the sphere. These points are illustrated in 
Fig. Q] where the bound-state and continuum densities are plotted for Al at metallic 
density and temperature k B T = 5 eV. Occupation numbers of bound states and con- 
tinuum partial-wave states inside the WS sphere are given along with bound-state 
eigenvalues in Table [T] 

The boundary conditions used here differ from those used by Sahoo et al. in 
Ref. 112511 . where the first derivative of the wave function is required to vanish at R vs . 
The differences in boundary conditions lead to major differences in the average- 
atom structure. For example, the model used in Ref. 12511 predicts that the M shell of 
metallic Al is partially occupied at temperatures k B T < 10 eV, whereas the present 
model predicts that the M shell is empty in this temperature range. Consequences of 
such differences are discussed later in Sec. |4] 
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3 Dynamic Structure Function 

As noted in the introduction, the theoretical model developed by Gregori et al. J^, 
with input from the average-atom model, is used here to evaluate the dynamic struc- 
ture function S(k,co). The ion-ion contribution Su(k,co) is evaluated in terms of 
Fourier-transforms of electron densities and the static ion-ion structure function 



Sa(k). Formulas for Su(k) are given in Ref. 13511 . Here, we use modified formulas 
that include options discussed in Ref. (H] for different electron and ion tempera- 
tures. The dominant effect of different electron and ion temperatures is to modify 
the relative size of elastic to inelastic contributions to S(k, co). The use of different 
temperatures is, therefore, a convenient tool for fitting experimental data, even in 
cases where equilibrium is reached. The inability to fit experimental data in some 
equilibrium cases without such an artifice is a weakness in the present scheme and 
is the subject of current research. The electron-electron contribution S ee (k, o) is ex- 
pressed in terms of the dielectric function e(k, co) of the free electrons which in turn 
is evaluated using the random-phase approximation (RPA). Plasmon resonances are 
present in S ee (k, co) at low momentum transfers k. Finally, bound-state contribu- 
tions to the dynamic structure function are evaluated using average-atom bound- 
state wave functions for the initial state. The final-state continuum wave function is 
described in two different ways: (1) approximating the final-state by a plane wave 
as in Ref. |25ll . and (2) using an average-atom final-state that approaches a plane 
wave asymptotically. There are dramatic differences between these two choices, 
especially at low momentum transfers. The more realistic average-atom choice au- 
tomatically includes ionic Coulomb-field effects. 



3.1 Ion-Ion Structure Function 

The contribution to the dynamic structure function from elastic scattering by elec- 
trons following the ion motion Su(k, co) is expressed in terms of the corresponding 
static ion-ion structure function Su(k) as: 

S u (k,co) = \f(k) + q (k)\ 2 Sn(k)8(co). (12) 

In the above, f(k) is the Fourier transform of the bound-state density and q(k) is 
the Fourier transform of the density of electrons that screen the ionic charge. In the 
average-atom approximation, the screening electrons are the continuum electrons 
inside the Wigner-Seitz sphere and 

r R ws , 

f(k)+q(k)=4n r 2 [n b {r)+n c ( r )}Mk r )dr, (13) 
Jo 

where jo (z) is a spherical Bessel functions of order 0. Note that f(Q) + q(Q) = Z in 
the average-atom model. In the applications discussed later, 8(co) is replaced by an 
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Fig. 2 Left panel: Sa(k) is shown for Be metal at electron temperature k B T e = 20 eV and ion- 
electron temperature ratios Tf/T e = (1, 0.5, 0.1) illustrated in solid, short-dashed and long dashed 
curves , respectively. The value k = 0.543 corresponds to an incident photon a>o = 2960 eV scattered 
at angle 40°. Right panel: Su(k, a>) in a.u. for Be metal at k B T e = 20 eV and Ti/T e = 0.1, where 
the function 8(a)) is replaced by a Gaussian of width 10 eV and C0[ = Wq — CO is the energy of the 
scattered x-ray. 



"instrumental" Gaussian, with full- width at half maximum =10 eV. This value is 
chosen because typical experiments in Be lfl3ll utilize a spectrometer with a 10-eV 
instrument width and use a CI Ly-a source at 2.96 keV. 

Approximate schemes to evaluate the static structure functions Sa(k) are dis- 
cussed, for example, in Ref. l36ll . Here, we follow Gregori et al. 0] and make use 
of formulas given by Arkhipov and Davletov l35ll that account for both quantum- 
mechanical and screening effects. The function Sa(k) in Ref. ll35tl is expressed in 
terms of the Fourier transform of the ion-ion interaction potential (r) through the 
relation: 

Su(k) = l--£ F &u(k). (14) 
K. B i 



Different Electron and Ion Temperatures 

In the average-atom model, T is the electron temperature T e which, in equilibrium, is 
equal to the ion temperature 7}. To allow for different electron and ion temperatures, 



the equations for Sn{k) given by Arkhipov and Davletov IB5I1 are modified following 
the prescription laid out by Gregori et al. [83. The electron temperature T e is replaced 
by an effective temperature T' e that accounts for degeneracy effects at temperatures 
lower than the Fermi temperature 7> . Similarly, the ion temperature 7] is replaced by 
an effective temperature T( that accounts for ion degeneracy effects at temperatures 
lower than the ion screened Debye temperature To- Explicit formulas for Su(k) are 
found in Ref. Ji| . The dramatic effect of different electron and ion temperatures on 
the static structure functions Sa(k) for metallic Be at T e = 20 eV are illustrated in 
the left panel of Fig. [2] This figure is similar to the upper-left panel of Fig. 1 in 
Ref. [8], which was obtained under similar condition. In the right panel of Fig. [2] 
contributions to Su(k, <X>) for Be at T e = 20 eV and T, ■ — 2 eV are shown. 
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3.2 Electron- Electron Structure Function 



The electron-electron structure function S ee (k, (o) is expressed in terms of the free- 
electron dielectric function e(k, (o) through the relation Eq. (42) in Ref. l27tl : 



See(k,0)) = - 



1 



1 -exp(-co/k B T) An 2 n, 



-Im 



1 



e(k,co) 



(15) 



We follow Gregori et al. y[ and evaluate the dielectric function in the random- 
phase approximation (RPA): 



, (16) 



where 



k 2 -2pkfJ. + 2(0 + iv k 2 + 2pkii-2<3-iv 



1 



(17) 



l+exv[(p 2 /2-n)/k B T] 

is the free-electron Fermi distribution. The parameter v in Eq. [16] is a positive in- 
finitesimal. The chemical potential ji is obtained from the average-atom model. It 
should be noted that the RPA dielectric function reduces to the Lindhard dielectric 
function J37, 38] in the limit k B T — > 0, where &(p) reduces to a unit step function 
that vanishes for p > k F . In the finite T case, we find 



djJ. 



i k 2 -2pkn + 2(0 + iv 
1 dn 



1 

2pl 
1 



In 



-l k z + 2pkfi-2co-iv 2pk 



In 



k 2 


+ 2pk + 


2co4 


iv 


k 2 


-2pk\ 


2co4 


iv 


k 2 


+ 2pk- 


2co- 


iv 


k 2 


-2pk- 


2(0- 


iv 



(18) 
(19) 



From the above expressions, it is clear that Re[e(k, (o)] is an even function of (O and 
that Im[e(£, ©)] is an odd function of (O. Therefore, we can limit ourselves to cases 
where (O is positive. For (O > 0, one obtains in the limit v — > 0+, 



lim s In 

v^O 



k 2 + 2pk + 2(0 + iv 
k 2 -2pk + 2(0 + iv 



^k 2 + 2pk-2(0-iv 



In 



k 2 -2pk-2(0-iv 



= ln 



k 2 +2pk + 2(0 



k 2 -2pk + 2(0 



In 



k 2 + 2pk-2(0 



k 2 -2pk-2(0 



- in only in the interval \k 2 - 2(o\ < 2pk < k 2 + 2(0. (20) 



It follows that 
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Fig. 3 Left panel: Real and Imaginary parts of the dielectric function e(k,a) are plotted along 
with — Im[l/e] for coherent scattering of a 2960-eV photon at 20° from Be metal at k B T = 18 eV. 
Right panel: The resulting electron-electron structure function S ee (k,a) (solid line) is shown to- 
gether with Sa(k,a>) (dashed line) which has an instrumental width of 10 eV. In this example 
k = 0.2757 a.u. and ffi>i = 2960 — CO (eV) is the energy of the scattered x-ray. 



Re[e(fc,©)] = l + ^3^ &{p)pdp 

k 2 + 2pk + 2co 



In 



k 2 -2pk + 2(0 



In 



k 2 + 2pk-2co 



k 2 -2pk-2(0 



and 



lm[e(k,co)} 



^(p)pdp 



2k B T 
k? 



■log 



l + esKp[Qi-a 2 /2)/k B T\ 



l+cxp[(fl-b 2 /2)/k B T] 



(21) 



(22) 



with a = \2co-k 2 \/2k and b = (2co + k 2 )/2k. From Eq. O and the fact that 
Im[l /e(k, co)] is an odd function of CO, it follows that 



S ee (k,eo) = exp(co/k B T)S ee (k,-co). 



(23) 



Therefore, in the absence of bound-state contributions, measuring the cross-section 
ratio at down-shifted and up-shifted plasmon peaks, provides an experimental 
method for determining the plasma temperature. 

The real and imaginary parts of e(k, co), along with — lm[l /e(k, £»)], are shown 
in the left panel of Fig. [3] for Be metal at temperature k B T = 18 eV. A plasmon 
peak, associated with coherent scattering of the incident x-ray by the electrons in 
the plasma, is seen near CO — 20 eV, where the real part of the dielectric function 
vanishes. The corresponding structure function S ee (k, co), illustrating up- and down- 
shifted plasmon features, is shown in the right-hand panel. For comparison pur- 
poses, the elastic-scattering structure function Su(k,co) is also shown in the right- 
hand panel. The functions displayed in this plot correspond to scattering of an inci- 
dent 2960-eV photon at an angle 20°. The classical plasma frequency for this case 
(n e = 2.05 x 10 23 cc- 1 ) is afo = 16.8 eV. 
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Fig. 4 Left panel: Real and Imaginary parts of the dielectric function e(k,a) are plotted along 
with — Im[l/e] for incoherent scattering of a 2960-eV photon at 90° from Be metal at k B T = 18 eV. 
Right panel: The electron-electron structure function S ee (k, co) (solid line) is shown together with 
Sa(k, Co) (dashed line) which is assumed to have an instrumental width of 10 eV. The corresponding 
momentum transfer is k = 1 . 1 23 a.u. and a>i = 2960 — ft) (eV) is the energy of the scattered x-ray. 



The amplitude and width of plasmon peaks is governed by the coherence param 



eter a = \/{X s k), defined in Eqs. (5-7) of Ref. 12311 . Here, X s is the shielding length 
given by 

/ k B TF l/2 (n/k B T) 
As -\j4nn e F_ 1/2 ^/k B T)' (24) 

where Fj(x) is a complete Fermi-Dirac integral, 



FJx) = — r / dy —, . (25) 

VV ; r(l + v)7o y l+exp{y-x) 

For the example shown in Fig. [5J the screening length A s = 1.440 and the cor- 
responding coherence parameter a = 2.520, so one expects and observes plasmon 
resonances. At values of a less than than one, coherent scattering by the plasma 
no longer occurs and the plasmon peak in the dielectric function disappears. An 
example of this behavior is illustrated in Fig. |4] where the dielectric function for 
scattering of a 2960-eV photon at 90° from Be metal at k B T = 18 eV is illustrated. 
The screening radius in this case remains unchanged A s = 1 .440, but the momentum 
transfer increases to k = 1.123 and the resulting coherence parameter is a = 0.6188. 
Therefore, as expected, all signs of plasmon resonances are absent in S ee (k, co). 



3.3 Bound-Free Structure Function 



The contribution to the dynamic structure function for Thomson scattering from 
ionic bound states Sb(k, co) is the sum over contributions from individual subshells 



X-ray Scattering from WDM 



11 



with quantum numbers (n,l): 

S b (k,CO) = £S„;(fe,t») 
nl 



Snl(k,CO) 



O n l 



21 



y f pd£2 p 



d 3 r\i/l(r)e ikr \l/„i m (r) 



(26) 
(27) 



e p =co+e„i 



where o nm is the fractional occupation number of subshell («,/), and k = ko — k\ 
and (0 = C0q — (0\ are the momentum and energy transfers, respectively, from the 
incident to the scattered photon as defined in Eq. ([T). In the above, y/p(r) * s an 
average-atom wave function that approaches a plane wave e ,p r asymptotically. We 
consider two cases below: firstly, the case where y p {r) is approximated by a plane 
wave and secondly, the case where y p (r) is an average-atom scattering function that 
approaches a plane wave asymptotically. 



Plane-wave final states 



Approximating the final state wave function by a plane wave Wp( r ) = e' pr , the 
bound-free structure function in Eq. ( f27l > can be rewritten as 



Snl(k,G>) 



O n l 



21 



pdQp 



d^re^WnUr) 



(28) 



where e p = (0 + £„/ and q = k~ p. Note that q is the momentum transferred to the 
ion. This expression may be simplified to 



S„i(k,C0) 



Onl 



p+k 



qdq\K„i(q)\ 



nk J\ p -k\ 

K m{q)= rdrji{qr)P n i{r). 
Jo 

S„i(k, (o) in Eq. ( f29b depends implicitly on ft) through the relation 



(29) 
(30) 



p= ^2(co + e„i 



Average-atom final states 

In the average-atom approach, the final state wave function consists of a plane wave 
plus an incoming spherical wave. (N.B. An outgoing spherical wave is associated 
with an incident electron. Time-reversal invariance, therefore, requires that a con- 
verging spherical wave be associated with an emerging electron.) The scattering 
wave function \j/p(r) in Eq. dZTT i may be written 
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Fig. 5 The beryllium S'-shell structure function Si s (k, Co) is shown for incident photon energy 
2960 eV and scattering angles 30° and 150°. The black curves show plane- wave results, the 
blue lines show the result obtained using an average-atom final-state wave function and the red 
lines show exact nonrelativistic Coulomb results. The dramatic suppression of average-atom and 
Coulomb structure functions at forward angles (the corresponding curves are multiplied by 10) is 
evident in the left panel. 



¥ P (r) = — L i h e- i5 h -P £h (r)Yr imi (p)Y limi (r), 



P i 



(31) 



where the continuum wave function is normalized to a phase-shifted sine wave 
asymptotically 

P sl (r) -+sm(pr-l7i/2 + 8i). (32) 
The factor e' k r in Eq. d27l i is expanded as 



«* r = 4;r£i^(fr)y^(*)y^(f). 

/2 m 2 



The bound-free structure function in Eq. ( l27b may then be expressed as 

c ° nl [ pdQp V 1 1 I 2 
"' 21 + 1 J (27r) 3 ^ | n,m| ' 

where 



with 



and 



hihM = -e l ^ p) J~ drP eh {r)P nl {r)j h {kr) 
{hm x \Y lm \hm 2 ) = JdQY,* mi (?)Y lm (r)Y hm (f). 



(33) 



(34) 



Jnlm = ^-Y, ih ~ h hlhM I (hm l \Y lm \l 2 m 2 )Y ll , ni (p)Y l * 2m2 (k), (35) 

P /,/ 2 m 1 m 2 



(36) 
(37) 
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Squaring /„/,„, integrating over angles of p, and summing over magnetic quantum 
numbers, one obtains 

Sni = — o„i2^Ai 1 ,i 2 \I h i h (p,k)\ (38) 
n hh 

where ^ 

A /l//2 = (2/ 1 + l)(2/ 2 +l)^ ( ;^ / 2 ^ . (39) 

In Fig. [5] several calculations of the structure function S\ s {k, Co) are compared 
for a photon of incident energy 2960 eV scattered at 30° (k = 0.41 1 a.u.) and 150° 
(k = 1.533 a.u.) from the K shell of beryllium metal at k B T = 20 eV. The results of 
calculations carried out using average-atom final states are smaller than those us- 
ing plane-wave final states by a factor of about 40 at 9 = 30° and 2.5 at = 150°. 
This suppression is a characteristic Coulomb field effect. Indeed, exact nonrelativis- 
tic Coulomb-field calculations of Thomson scattering l39ll . with nuclear charge ad- 
justed to align the Coulomb and average-atom thresholds, illustrated by the curves 
labeled "Coulomb" in Fig. [5] show a similar suppression. 

The situation is quite different at much higher momentum transfer. This fact is 
illustrated in the left panel of Fig. [6] where we consider scattering of an 18-keV 
photon at 130° (k = 8.750 a.u.) from Be metal at k B T = 12 eV. The solid black curve 
shows the bound-state contribution to the S(k, co) evaluated using an average-atom 
final state, while the dashed green line shows the contribution evaluated using a 
plane-wave final state. In this high-momentum transfer case, the two calculations 
differ qualitatively near threshold and the peak energies differ by a few percent, but 
the average-atom and plane-wave results are otherwise quite similar. 

In the right panel of Fig. [6] we consider a more complex example, scattering of a 
9300-eV x-ray at 130° (k = 4.521 a.u.) from tin metal at k B T = 10 eV. Contributions 
from individual subshells are shown together with their sum Si,(k, co). 



4 Applications 

The scheme developed in the previous sections is applied to several examples of 
current or possible future interest. The example of hydrogen at electron density 
10 24 cc _1 and temperature 50 eV is considered first. In this example, the H ion is 
completely stripped and inelastic scattering is determined completely by S ee (k, Co). 
The plasmon resonances present in the hydrogen spectrum at forward scattering an- 
gles fade rapidly with increasing angle. The case of beryllium metal at 18 eV, which 
has been well studied both theoretically and experimentally, is considered next and 
good agreement is found between average-atom predictions and experimental data. 
For beryllium, the bound-state contribution is small and beyond the range of avail- 
able experimental data. As a third example, calculations of S(k, co) for scattering of 
9.3-keV and 3.1-keV x-rays from aluminum metal at k B T = 5 eV are evaluated and 
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Fig. 6 Left panel: The dynamic structure function for the Is state of Be metal at k B T = 12 eV 
is shown for scattering of an incident x-ray with energy coq = 18-keV x-ray at 130°. The solid 
black and dashed green lines are results calculated using average-atom and plane-wave final states, 
respectively. The dashed red line is the elastic scattering contribution Sa(k, ft)) with instrumental 
width 10 eV. Right panel: The total bound-free dynamic structure function Si,(k, ft)) with average- 
atom final states is shown for tin (Z = 50) at metallic density and k B T = 10 eV is shown, along 
with contributions from individual subshells. The incident x-ray energy is ff^ = 9.3 keV and the 
scattering angle is 130°. In the above, a>\ = ft)o — ft) is the energy of the scattered x-ray. 



compared with the average-atom predictions by Sahoo et al. 112511 . Finally, Thomson 
scattering of 9.3-keV x-rays from titanium and tin at 30° and 130° are considered 
to illustrate cases where bound-state contributions are the dominant features of the 
scattered x-ray spectrum. In Table [2] we list some important average-atom proper- 
ties for the elements considered in the following subsections. 



4.1 Hydrogen 

In the average-atom model, a density p = 1.931 g/cc is required at k B T = 50 eV 
to achieve free-electron density n e — 10 24 cc _1 . Physical properties of the plasma 
under these conditions are listed in the second column of Table [2] The single elec- 
tron in hydrogen is completely stripped, leaving only the continuum inside the WS 
sphere. The continuum density n c (r) merges into the uniform free-electron den- 



Table 2 Average-atom parameters for the examples of hydrogen, beryllium, aluminum, titanium, 
and tin presented in Sec. [4] 





H 


Be 


Al 


Ti 


Sn 


k B T (eV) 


50 


18 


5 


10 


10 


P (g/cc) 


1.931 


1.635 


2.700 


4.540 


7.310 


R ws («o) 


1.118 


2.452 


2.991 


3.044 


3.515 


H (a.u.) 


1.091 


-0.531 


0.241 


-0.044 


-0.071 


Abound 





1.966 


10.000 


17.341 


45.622 




1 


2.036 


3.000 


4.659 


4.378 


2/ 


0.867 


1.647 


2.146 


2.322 


3.374 


n e (10 23 cc _1 ) 


10 


1.800 


1.292 


1.326 


1.251 
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r (a.u.) «! (keV) 

Fig. 7 Left panel: The continuum density n c (r) inside the WS sphere for hydrogen at k B T = 50 eV 
and density 1.931 g/cc is seen to converge smoothly to the free-electron density n e = 10 24 cc~' 
beyond the WS radius R ws . Right panel: The dynamic structure function S(k, ft)) for H at k B T — 
50 eV and n e = 10 24 cc~' is plotted for scattering of a 5-keV photon at angles 20°, 30° and 40°, 
and ft>i = 5 — ft) (keV) is the energy of the scattered x-ray. 



sity n e outside the sphere. The total number of electrons inside the WS sphere 
N c = 4;r/ ws r 2 n c {r)dr = 1, however, the number of free electrons per ion is 
Zf = 0.867. These points are illustrated in the left panel of Fig, [7] 

Since there are no bound electrons, only S ee (k,co) contributes to the inelastic 
scattering. The dynamic structure functions for scattering of a 5-keV x-ray at angles 
20°, 30° and 40° are shown in the right panel of Fig. UJ The prominent plasmon 
peaks at 9 = 20° wash out at as the scattering angle (and momentum transfer k) 
is increased. For this particular case, the screening length X s = 1.071 a.u. differs 
only slightly from the WS radius /? ws = 1.118 a.u.. The values of the coherence 
parameter a for = (20°, 30°, 40°) are a = (2.005, 1.345, 1.018), respectively. 
The resonant features in Fig. |7jare seen to be distinct for a > 1 but disappear as 
a approaches 1. It should be noted that the (unperturbed) plasma frequency for 
hydrogen at n e = 10 24 cc~' is co p i = 37.1 eV. 



4.2 Beryllium 

In the left panel of Fig. [8] the dynamic structure function for scattering of a 2963- 
eV photon at 40° from beryllium (density = 1.635 g/cc) at k B T = 18 eV is shown. 
The L-shell electrons are completely stripped under these conditions but the K shell 
remains 97% occupied. The ion temperature, which governs the amplitude of the 
elastic peak, is chosen to be 2. 1 eV in this example. For the case at hand, the coher- 
ence parameter is a = 1.21 > 1, so one expects and observes plasmon peaks in the 
scattering intensity profile. In the average-atom approximation, the threshold energy 
for bound-state contributions is the average-atom eigenvalue ei s = —86.8 eV. (As 
mentioned earlier, average-atom eigenvalues are inaccurate approximations to re- 
moval energies. The average-atom threshold \e\ s \ is 20% smaller than the measured 
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2840 2880 2920 2960 3000 2880 2920 2960 3000 

(Bj (eV) cOj (eV) 

Fig. 8 Left panel: Structure function S(k, co) for scattering of a 2963-eV photon at 40° from beryl- 
lium at n e = 1 .8 x 10 23 cc~' and electron temperature 18 eV. The bound-state contribution is scaled 
up by a factor of 50. Right panel: Intensities (arbitrary units) for scattering of a CI Ly-a x-ray from 
beryllium at 40° l4lll ; measured data (black solid line), source function (blue dashed line), and the 
average-atom fit (red solid line). In the above, oj] is the energy of the scattered x-ray. 



^T-shell threshold 1 1 1 eV in beryllium metal J4(J.) In the average-atom approxima- 
tion, bound-state contributions to S(k, co) from ^T-shell electrons have a threshold at 
COi = 2963 - 86.8 = 2876.2 eV. The K-sheU contribution multiplied by 50 is shown 
in the left panel. The average-atom parameters used in this calculation are listed in 
column three of Table 12 

To validate the present average-atom model against experimental data, a Be ex- 
periment done at the Omega laser facility that used a CI Ly-a source to scatter from 
nearly solid Be at an angle of 40° is used. An electron temperature of 18 eV, ion tem- 
perature of 2.1 eV, and density of 1.635 g/cc used in the average-atom model give s 
an electron density of 1.8 x 10 23 cc~', in agreement with the analysis in Ref. 114 ill . 
The right panel of Fig.[8]shows the experimental source function from the CI Ly-a 
line as a blue dashed line. Because of satellite structure in the source we approx- 
imate the source by three lines: a CI Ly-a line at 2963 eV with amplitude 1 and 
two satellites at 2934 and 2946 eV with relative amplitudes of 0.075 and 0.037, 
respectively. Doing the Thomson scattering calculation using the three weighted 
lines, we calculate the scattering amplitude for Thomson scattering (red solid line) 
and compare against the experimental data (black solid line) here. We observe ex- 
cellent agreement within the experimental noise. Contributions from the bound Is 
electrons, which have a threshold at 2876 e V, are beyond the range of the data shown 
in the right panel. 



4.3 Aluminum 

As mentioned in the introduction, the present average-atom calculations disagree 
in various ways with those in Ref. Ir25ll . where Thomson scattering from aluminum 
metal at several temperatures was considered. To clarify the differences, we compare 
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results of the present calculations for the case of aluminum at k B T = 5 eV with the 
corresponding results presented in Ref. l25ll . The dimensionless product S(k, ft)) (0p\ 
is shown for scattering of 9300-eV and 3100-eV x-rays in the left and right pan- 
els of Fig. [9] respectively. The upper and lower panels of Fig. [9] show results for 
scattering angles of 30° and 130° respectively. The plots in Fig.|9]can be compared 
directly with those in Fig. 6 and Fig. 7 of 12511 . where similar plots for 9300-eV 
and 3100-eV x-ray energies can be found. The size and shape of the free-electron 
contributions shown in Fig. [9] are in agreement with those shown in the correspond- 
ing figures in [25], however, the L-shell contributions to S(k, ft)) in i25ll are larger 
by a factor of approximately 5 than those shown in Fig. [9] Moreover, contributions 
from the M shell, which dominate the spectrum just below the elastic peak in l25ll 
are completely absent in Fig. [9] Owing to the differences in bound-free contribu- 
tions Sb(k, ft)), the present prediction for the aluminum dynamic structure function 
S(k, ft)) differs in both size and shape from that given in 12511 . 

Differences in boundary conditions between the two average-atom models ac- 
count for the presence or absence of bound M-shell electrons as discussed in Sec. [2] 
Furthermore, difference in the relative size of the L-shell contributions at forward 
angles may owe in part to the use of plane-wave final states in 1E5I1 . However, the 
substantial difference in size of the L-shell contributions at backward angles remains 
a mystery. The average-atom parameters for aluminum used in the present calcula- 
tion are given in the fourth column of Table [2] 



4.4 Titanium & Tin 

The average-atom model predicts that a titanium metal at 10 eV is in an Ar-like 
configuration with completely filled K and L shells together with 1 .97 3s electrons 
5.36 3p electrons in the M shell. This is an interesting case where contributions to 
St,{k,(o) from the L and M shells are substantial. In the left panels of Fig. [TOl the 
dynamic structure function S(k, ft)) is shown in the solid black line for an incident 
9.3-keV x-ray scattered at 30° and 130°. Contributions from Sb(k, ft)) shown in the 
long-dashed green lines dominate the free-electron contributions S ee (k,(o) shown 
in the short-dashed red lines. As seen in the lower-left panel, sharp features associ- 



Table 3 Comparison of the theoretical average-atom thresholds (—£„/) and experimental thresh- 
olds from the NIST database |40| | for contributions to S/, (k, co) from individual subshells in titanium 
and tin. 



Titanium Tin 



nl 


theory 


expt. 


A% 


nl 


theory 


expt. 


A% 


2s 


512.7 


563.7 


9 


3d 


477.8 


488.2 


2 


2p 


426.8 


457.5 


1 


As 


118.9 


136.5 


13 


is 


45.4 


60.3 


25 


*P 


83.1 


86.6 


4 


3p 


22.8 


34.6 


34 


4cl 


22.4 


23.9 


6 










5s 


2.1 
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Fig. 9 S(k, a) Opi is shown for scattering of 9300-eV and 3100-eV x-rays at 30° and 130° from Al 
metal at k B T = 5 eV. The short-dashed red curves are contributions S ee (k, CO) from free electrons, 
the long-dashed green curves give contributions Sf,(k, Co) from bound 2s and 2p electrons and the 
solid black curves show the total dynamic structure function. Instrumental widths of 0.1 %C0q were 
used. For Al metal at k B T = 5 eV, the plasma frequency is dOp] = 13.35 eV. In the above figure, tt>i 
and cat) are energies of the scattered and incident x-rays, respectively, and CO = COq — CO\ . 



ated with excitation of the M subshells (3p and 3s) show up just below the elastic 
peak at a>\ = (Oq, whereas features associated with excitation of the L subshells {2p 
and 2s) show up 400 to 500 eV below the elastic peak. The bound-state thresholds 
for titanium are compared with measured thresholds from the National Institute for 
Standards & Technology (NIST) database S in Tabled 

Metallic tin at 10 eV has a Pd configuration with 45.6 bound electrons. The 
resulting dynamic structure function S(k, co) is shown in the right panels of Fig. [10] 
for scattering angles of 30° and 130°. The situation is similar to that for titanium; 
bound-free contributions dominate those from free-electrons. The irregularities in 
the bound-state contributions are associated with the onset of contributions from 
individual subshells. For scattering at 130°, contributions from individual subshells 
are shown in the right panel of Fig. [6] Average-atom thresholds for tin are compared 
with values from the NIST database H in Tabled 
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Fig. 10 Left panels: S(k, ft)) for scattering of a 9.3-keV x-ray at = 30° and 130° from Ti metal 
at kg = 10 eV is shown in the solid black line. The free-electron contribution S ee (k, ft)) is shown in 
the short-dashed red line while the contribution from bound states Si, (k, ft)) is shown in the long- 
dashed green line. Right panel: Same as left panel, except scattering from Sn metal at k B = 10 eV. 
Contributions from individual subshells of Sn are shown in the right panel of Fig.flO] In the figure 
above, ft>i = 9300 — ft) (eV) is the energy of the scattered x-ray. 



5 Summary 

A scheme for analysis of Thomson scattering from plasmas based on the average- 
atom model is presented. Given the plasma composition (Z,A), density p and tem- 
perature T, the model gives, in addition to the equation of state of the plasma, all 
parameters needed for a complete description of the Thomson scattering process. In 
particular, the average-atom code predicts wave functions for bound and continuum 
electrons, densities of bound, screening, and free electrons, and the chemical poten 



tial. Predictions of the present average-atom model disagree with those in Ref. 125 ] 



where a similar model with different boundary conditions was used. In the average- 



atom model used in Ref. 12511 . 3d (M shell) electrons are bound in metallic Al for 
temperatures between 2 and 10 eV, leading to substantial bound-state contributions 
to the dynamic structure function. In the present model, by contrast, the M shell of 
metallic aluminum is vacant in the temperature range k B T < 10 eV and the corre- 
sponding bound-state features are absent. 
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Elastic scattering from bound and screening electrons is treated here following 
the model proposed by Gergori et al. J3] which makes use of formulas for the static 
ion-ion structure function Su{k) given by Arkhipov and Davletov ll35ll . Modifica- 
tions suggested in Ref. JU to account for different electron and ion temperatures are 
also included. Specifically, in the applications considered here, the amplitude of the 
elastic peak is adjusted artificially by choosing an ion temperature that is different 
from the electron temperature, even in cases where equilibrium is expected. Such an 
adjustment was used to fit the experimental data for beryllium shown in Fig.[H] The 
dynamic structure function for scattering from free electrons depends sensitively on 
the free-electron dielectric function e(k, co). Again, we follow the model proposed 
in Ref. ^ and evaluate the dielectric function in the random-phase approximation. 
The RPA dielectric function includes features such as plasmon resonant peaks that 
show up in experimental intensity profiles and can be used in connection with the 
principle of detailed balance to determine electron temperatures. Bound-state fea- 
tures are easily included in the present scheme, inasmuch as the average-atom model 
provides bound-state and continuum wave functions. Ionic Coulomb-field effects 
are features of calculations carried out using average-atom wave functions rather 
than plane waves to describe the final state electrons. In conclusion, the average- 
atom model provides a simple and consistent point of departure for the theoretical 
analysis of Thomson scattering from plasmas. 
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